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I.  INTRODUCTION 


Mgh  range  resolution  has  maity  advantages  in  radar.  Apart  from  providing  the  ability 
to  resolve  closed  spaced  targets  in  range,  it  improves  the  range  accuracy?,  reduces  the  amount 
ofclutter  within  the  resolution  cell,  reduces  multipath,  provides  high  resolution  range  profiles 
and  aids  in  target  classification.  Typical  waveforms  for  achieving  high  resolution  are  ultra 
wideband  (UWB),  conventional  pulse  compression  and  step  frequency  waveforms.  In  all  three 
waveforms,  high  range  resolution  is  obtained  by  increased  bandwidth.  In  the  first  two 
waveforms  large  bandwidth  is  instantaneous,  whereas  in  the  step  frequency  waveform  large 
bandwidth  is  achieved  sequentially  over  many  pulses  while  keeping  the  instantaneous 
bandwidth  low.  This  results  in  lowering  of  A/D  sampling  rate  for  the  step  frequency 
waveform.  This  is  a  significant  advantage  as  limitation  of  A/D  technology  is  a  bottleneck  in 
the  implementation  of  high  resolution  radar. 

The  use  of  a  step  frequency  waveform  for  detailed  RCS  measurements  in  anechoic 
chambers  or  in  opOT  ranges  is  well  known.  In  these  applications,  the  step  frequency  waveform 
along  with  ISAR  processing  provides  high  resolution  mapping  of  rotating  targets.  The  use  of 
the  step  frequency  waveform  for  operational  applications  needs  to  be  explored.  The  purpose 
of  this  report  is  to  design  the  step  frequency  waveform  for  detection  of  small  targets 
embedded  in  clutter.  Any  solution  to  this  challenging  radar  problem  would  include  large 
improvements  in  signal-to-noise  (S/N)  and  signal-to-clutter  (S/C)  ratios.  S/N  ratio  is 
improved  by  having  large  average  transmit  power  or  integrating  large  number  of  pulses  or 
some  combination  of  both.  S/C  ratio  is  improved  by  reduction  of  clutter  which  is  achieved  in 
two  steps,  that  is,  first  by  limiting  the  amoimt  of  clutter  entering  the  radar  receiver,  and 
secondly  by  canceling  the  clutter  in  doppler  processing.  The  amount  ofclutter  entering  the 
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radar  receiver  is  achieved  by  decreasing  the  effective  pulse  vwdth  which  is  normally  done  with 
a  conventional  pulse  compression  waveform  to  keep  adequate  S/N  ratio. 

It  should  be  noted  that  conventional  constant  frequency  medium  or  high  PRF 
waveforms  would  satisfy  the  above  requirements  for  large  improvement  in  S/N  and  S/C 
ratios.  However,  the  constant  frequency  waveforms  would  require  high  speed  A/D’s  to 
achieve  small  range  cell  Hi^  PRF  waveforms  will  be  highly  ambiguous  in  range  and  the 
range  ambiguities  will  ftuther  increase  the  clutter  problem  as  clutter  from  various  range  zones 
due  to  different  pulses  arrive  simultaneously  to  compete  with  the  target  signal.  To 
accommodate  the  large  folded  clutter,  A/D  would  require  a  large  dynamic  range  and  high 
sampling  rate,  a  combination  difficult  to  achieve  in  practice.  Thus,  requirements  for  systems 
with  fine  range  resolution  may  not  be  met  with  the  existing  A/D  technology.  Step  frequency 
can  alleviate  the  A/D  problem  by  easing  the  requirements  both  on  the  sampling  rate  and  the 
dynamic  range. 

Range  resolution  in  step  frequency  is  controlled  by  frequency  spread  over  the  pulse 
burst  (which  is  NAf,  where  N  is  the  number  of  pulses  in  the  burst  and  Af  is  the  frequency 
step)  but  the  sampling  rate  depends  upon  the  pulse  width.  Thus,  in  the  step  frequency  radar 
sampling  rate  is  made  independent  of  the  range  resolution  which  allows  lowering  its  value 
without  sacrificing  range  resolution.  The  dynamic  range  in  step  frequency  radar  is  reduced 
by  limiting  the  folded  clutter  from  ambiguous  ranges.  Unlike  constant  frequency  radars, 
multiple  time  around  clutter  in  step  frequency  radar  does  not  come  from  the  entire  range 
extent  but  is  limited  by  the  number  of  frequency  lines  which  will  pass  through  the  IF  filter. 
The  number  of  frequency  lines  and  thus  the  number  of  ambiguous  range  zones  is  1/tAf  (where 
T  is  the  pulse  width).  With  judicious  choice  of  waveform  parameters  it  is  possible  to  limit  the 


2 


dynamic  range  of  the  clutter  as  well  as  the  sampling  rate  of  AID  and  still  meet  the  high 
resolution  requirements. 

The  design  of  constant  fiequency  waveforms  is  relatively  simple  and  well  documented. 
However,  the  deagn  of  step  frequency  waveforms  is  complex  due  to  the  conflicting  effect  of 
waveform  parameters  on  various  desired  performance  goals.  Thus,  in  this  report  effort  is 
made  to  ^stematize  the  design  process  for  the  high  PRF  step  frequency  waveform  The  goals 
in  the  deagn  of  a  high  PRF  waveform  include  keeping  the  target  in  clutter  free  area,  keeping 
the  target  doppler  unambiguous  and  meeting  other  radar  performance  specifications. 

Also,  use  is  made  of  the  ambiguity  function  to  determine  the  characteristics  and 
special  advantages  of  the  step  frequency  waveform.  Finally,  ambiguity  diagrams  are  sketched 
for  a  specific  waveform  of  interest. 
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II.  INTRODUCTION  TO  THE  STEP  FREQUENCY  RADAR 


This  ch^ter  gives  a  brief  introduction  to  the  concepts  and  implementation  of  the  step 
frequency  radar  before  the  discussion  of  the  waveform  design,  which  is  done  in  the  next 
chapter. 

A.  STEP  FREQUENCY  WAVEFORM 

The  step  frequency  waveform  consists  of  a  series  of  N  pulses  each  with  a  pulse  width 
of  T,  and  whose  frequency  is  iira-eased  from  pulse  to  pulse  in  steps  of  Af,  as  shown  in  Figure 
1.  The  frequency  of  the  w***  transmitted  pulse  is  given  by 

=/o+(«-l)A/,  (2.1) 

where.4  is  the  nominal  carrier  frequency,  A/is  the  frequency  step  size  and  w  =  1, 2,  3, ...,  N. 


/q  “  Nominal  carrier  frequency 
A/“  Frequency  step  size 


PRI 


Figure  1.  Step  Frequency  Waveform. 
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Although  the  instantaneous  bandwidth  of  a  step  frequency  waveform  is  approximately 
1/t,  a  much  larger  effective  bandwidth  can  be  achieved.  This  effective  bandwidth  is  given  by 

=  W-  (2.2) 


The  range  resolution  for  this  waveform  is  not  the  conventional  value  of  ct/2.  The 
synthetic  range  resolution  obtained  by  coherently  processing  a  series  of  N  returns  from  the 
target,  can  be  much  smaller  and  is  given  by 


M  =  —  =  — , 
2B^  2AA/ 


(2.3) 


where  c  is  the  speed  of  light  (3xl0*m/sec).  The  fact  that  in  the  step  frequency  radar, 
resolution  does  not  depend  on  the  instantaneous  bandwidth  and  the  fact  that  resolution  can 
be  increased  arbitrarily  by  increasing  NAf,  are  sigmficant  advantages.  There  is  a  constraint 
on  selection  of  Af  (i.e.,  Af<  1/t),  however,  N  can  be  increased  to  achieve  high  range 
resolution. 

B.  SYSTEM  DESCRIPnON  AND  IMPLEMENTATION 

The  implementation  of  the  step  frequency  radar  is  similar  to  the  coherent  Pulse- 
Doppler  radar.  The  system  implementation  can  be  observed  in  Figure  2.  The  core  of  the 
system  is  a  coherent  step  frequency  S5mthesizer  with  an  output  frequency  called^^  which 
is  stepped  from  pulse  to  pulse  by  a  fixed  step  size,  A/  Therefore,  we  can  write 

fsyn  =  («-l)A/.  (2.4) 
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The  frequency  of  a  local  coherent  oscillator  (COHO),XoAo>  is  first  mixed  with  the  output 
frequency  of  the  step  frequency  synthesizer,^ .  Then  the  sum  frequency  is  converted  to  the 
final  transmitted  frequency  by  mixing  it  with  the  frequenqr  of  the  stable  local  oscillator 
(STALO),  .  The  sum  frequency  of  this  second  mixer  is  then  pulse  modulated  and 
amplified  before  the  transmission  occurs.  Each  transmitted  pulse  will  have  a  carrier  frequency 
comprised  of  three  different  conqjonents:  the  fixed  IF  frequency,  fixed  RF  frequency, 

Xato  and  the  variable  frequency  of  the  stepped  frequency  synthesizer,^  .  Thus,  the  carrier 
frequency  of  the  w*  pulse  is  given  by  the  following  equation; 

fn  =  fmlo^fcoho-^fsyn  =  /o+  («“  O  A/.  (2.5) 

In  the  receiver  side,  the  return  signal  is  first  down  converted  by  mixing  it  with_;^^/„ . 
The  resulting  agnal  is  fiuther  down  converted  to  intermediate  frequency  (IF)  by  mixing  with 
.  The  output  signal  obtained  after  the  second  mixer  is  the  IF  signal.  The  signal  is  then 
divided  into  two  different  channels,  the  in-phase  channel  (I),  and  the  quadrature  channel  (Q). 
Both,  the  in-phase  and  the  quadrature  signals  are  mixed  with  the  output  signal  of  the 
coherent  oscillator.  The  first  is  mixed  directly  and  the  second  is  mixed  with  the  COHO  signal 
after  this  last  one  has  been  shifted  90  degrees.  The  I  and  Q  components,  both  in  the  video 
frequency  range,  are  sampled  by  an  A/D  converter  at  a  rate  equal  to  the  inverse  of  the  pulse 
width. 

Each  complex  sample  is  usually  called  rcmge  bin,  thus  range  bin  may  be  defined  as  a 
memory  location  for  ten:q)orary  storage  of  successive  pairs  of  numbers  representing  the  I  and 
Q  samples  of  the  radar  return  received  at  a  given  point  in  the  interpulse  period.  A  separate 
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bin  therefore  must  be  provided  for  each  sampling  interval  (range  gate).  To  the  extent  that 
range  is  unambiguous,  the  numbers  stored  in  any  one  bin  represent  successive  returns  from 
a  single  range  increment,  hence  the  name  range  bin.  Because  of  the  correspondence  of  the 
range  bins  to  the  sampling  intervals  (when  A/D  conversion  follows  I  and  Q  detection),  range 
bin  has  come  to  be  used  synonymously  with  sampling  interval  as  well  as  range  gate. 

Sampled  returns  from  incoming  pulses  are  stored  im  memoiy  until  all  the  pulses  within 
the  same  burst  have  been  received  and  can  be  processed  as  shown  in  Figure  3.  Figure  4 
describes  the  organization  of  sample  storage  in  memory.  Complex  samples  for  each  range  bin 
are  transformed  by  taking  the  FFT  of  samples  in  each  range  bin,  resulting  in  a  high  resolution 
range  profile  (HRR  profile). 

One  important  advantage  of  the  step  frequency  radar  compared  with  other  radars 
which  also  use  wideband  waveforms,  is  that  in  the  former  the  narrow  instantaneous  bandwidth 
does  not  require  a  high  analog  to  digital  sampling  rate  which  can  be  a  major  limitation  in  the 
system  design.  On  the  other  hand,  the  step  frequency  radar  may  require  more  complex  signal 
processing. 


Figure  3.  Elementary  processing. 
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Figure  4.  N  samples  from  a  target  are  collected  and  stored  in  memory  for 
processing. 
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C.  WAVEFORM  PROCESSING 

The  radar  return  signal  at  the  output  of  the  coherent  detector  can  be  written  as 

sin)  =  Iin)+jQin)  =  ^^expCytf)^)  ,  (2.6) 

where 

A„  =  the  magnitude  of  the  m***  pulse; 

(f>„  =  the  phase  of  the  pulse. 

A„  represents  an  amplitude  factor  which  depends  on  the  transnutted  power,  the  size 
of  the  target,  and  other  factors  in  radar  range  equation. 

The  phase  of  the  /impulse  can  be  written  as 


4>„  =  (2.7) 

where/„  is  the  carrier  frequency  of  the  n*  pulse  and  is  the  round  trip  time  from  the  target 
measured  by  the  n  pulse.  This  time  is  equal  to  '2RJc ,  with  being  the  range  between  the 
target  and  the  radar.  For  a  target  with  constant  velocity  toward  the  radar,  the  range  K„  is 
given  by 

=  ^o-("-l)vr,  (2.g) 

where  Rp  =  initial  target  range; 

T  =  pulse  repetition  interval; 

V  =  relative  radial  velocity  between  target  and  radar. 
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The  overall  expression  for  the  baseband  return  signal  can  now  be  determined  by 
substituting  the  instantaneous  frequency  from  Equation  2.1,  and  the  target  ran^e  from 
Equation  2.8  in  Equation  2.7,  and  then  substituting  the  resulting  expression  in  Equation  2.6. 
The  resulting  equation  is  as  foUows: 


s(n)  =  A„  exp 


.  471 


=  exp  y 


.  47: 


fnK 


(2.9) 


Equation  2.9  can  be  rewritten  as 


s(n)  =  ^„exp 


+  27t^— («-l)r  +  27C-^(M-l)r+  2it  {n-\)T 


(2.10) 

The  first  term  represents  a  constant  phase  shift  which  is  not  of  any  practical 
significance.  The  second  term  is  the  multiphcation  of  the  rate  of  change  of  frequency  Af/T 
with  the  roimd  trip  time  ZR/c.  This  term  represents  a  shift  in  frequency  during  the  round  trip 
time  2R/c.  Thus  the  range  (or  the  round  trip  time)  is  converted  into  a  frequency  shift  (which 
is  analogous  to  conversion  of  range  to  frequency  in  linear  frequency  modulated  CW  radar). 
Therefore,  it  is  possible  to  resolve  and  measure  the  range  by  resolving  the  frequency,  which 
can  be  done  by  taking  the  DFT  of  the  recdved  agnal  from  N  frequency  stepped  pulses.  Since 
the  range  is  measured  by  taking  the  DFT,  the  range  measurement  will  have  the  same 
limitations  as  the  frequency  measurement  by  DFT.  Thus,  the  range  resolution  AR  and 
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unambiguous  range  are  dependent  on  the  frequency  resolution  and  the  mavimnm 
unambiguous  frequency  measurement  by  DFT,  respectively.  This  range  resolution  oic/2NAf 
is  equivalent  to  the  resolution  obtained  by  conventional  pulse  compression  with  a  compression 
ratio  ofNvAf.  The  maximum  value  of  the  compression  ratio  is  limited  to  N  tor  practical  cases 
(as  tAf  is  constrained  to  be  equal  to  or  less  than  unity).  For  detection  of  moving  targets, 
range  resolution  can  be  traded  off  for  creation  of  clear  space  within  the  range  window  R„.  The 
third  term  in  Equation  2.10  represents  the  doppler  frequency  shift  due  to  target  motion  and 
it  adds  to  the  frequency  shift  of  the  second  term.  The  range  resolution  process  unwillingly 
treats  the  doppler  frequency  as  a  frequency  shift  due  to  range  and  thus  results  in  shifting  of 
target  range  from  its  tme  range.  The  fourth  term  in  Equation  2. 10  is  due  to  the  interaction  of 
the  changing  frequency  of  the  step  frequency  waveform  with  the  target  motion.  The  doppler 
shift  changes  with  each  pulse  (even  for  constant  velocity  targets)  because  of  the  change  in 
pulse  frequency.  The  doppler  shift  due  to  the  constant  frequency  component  yj,  has  already 
been  taken  into  account  in  the  third  term.  The  fourth  term  gives  the  doppler  shift  2vnAf/c  due 
to  the  frequency  step  of  the  w*  pulse.  Thus,  the  return  from  a  moving  target  due  to  N 
frequencies  will  contain  N  frequency  components  in  the  data  domain  instead  of  ideally 
containing  one  component.  This  spread  in  frequency  of  2NvAf/c  when  processed  by  DFT  will 
lead  to  a  spread  in  range  by  vAT,  which  in  terms  of  range  bins  is  given  by 


III.  STEP  FREQUENCY  WAVEFORM  DESIGN 


Radar  waveforms  can  be  classified  as  low  PRF,  high  PRF  or  medium  PRF.  A  low  PRF 
is,  by  definition,  one  for  which  the  maximum  range  the  radar  is  designed  to  handle  lies  in  the 
first  range  zone  -  the  zone  from  which  first-time-around  echoes  are  received.  Range  for  low 
PRF  is  unambiguous  and  doppler  fi-equency  is  ambiguous.  A  high  PRF  is  one  for  which  the 
observed  doppler  frequencies  of  all  significant  targets  are  unambiguous  and  range  is 
ambiguous.  Finally,  a  medium  PRF  is  one  for  which  neither  of  these  conditions  is  satisfied. 
Both  range  and  doppler  fi-equency  are  ambiguous.  In  this  report  the  interest  is  in  the  high  PRF 
waveform. 

The  purpose  of  this  chapter  is  to  develop  a  design  procedure  for  detection  of  small 
targets  with  a  surface  (land  or  sea)  based  step  fi-equency  radar  which  employs  a  high  PRF 
waveform.  The  disadvantage  of  high  PRF  is  that  the  range  ambiguity  has  to  be  resolved. 
However,  several  advantages  result  fi-om  high  PRF.  The  unambiguous  doppler  leads  to 
advantages  such  as  elimination  of  clutter  and  absence  or  minimization  of  the  blind  doppler 
problem.  It  also  improves  the  signal-to-noise  (S/N)  ratio  by  integration  of  larger  number  of 
pulses.  For  a  constant  fi-equency  high  PRF  waveform,  multiple  time  around  clutter  (fi-om 
ambiguous  regions  due  to  different  pulses)  adds  up  and  thus  decreases  the  signal-to-clutter 
(S/C)  ratio  and  increases  the  dynamic  range  of  the  receiver  and  associated  A/D.  Though  the 
S/C  ratio  can  be  improved  by  doppler  filtering,  receiver  and  A/D  will  still  require  a  large 
dynamic  range  as  a  large  magnitude  of  folded  clutter  has  to  be  accommodated  in  the  receiver 
and  A/D.  The  step  frequency  waveform  reduces  the  folded  clutter  and  thus  the  dynamic  range 
by  decreasing  the  number  of  range  ambiguities.  Range  ambiguities  in  the  step  frequency 
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waveform  depend  upon  the  number  of  pulses  which  will  pass  through  the  IF  filter 
simultaneously.  This  number  is  ^ven  by  1/(tA/).  By  proper  choice  of  the  product  xA/  range 
ambiguities  and  thus  clutter  and  dynamic  range  can  be  reduced.  However,  this  advantage 
comes  about  at  the  cost  of  limiting  the  range  coverage.  This  feature  of  the  step  frequency 
waveform  provides  a  good  tradeoff  as  high  resolution  radars  are  often  limited  by  A/D 
technology.  The  waveform  design  involves  the  selection  of  PRF  (/),  number  of  pulses  (N), 
frequency  step  (Af)  and  pulse  width  (x).  It  is  assumed  that  the  following  radar  performance 
related  parameters  are  available:  range  resolution  (AR),  nominal  carrier  frequency  time- 
on-target  (tot),  minimum  and  maximum  radial  velocity  (v^  and  and  the  maximum 
number  of  pulses  (N^.  The  design  of  constant  frequency  waveforms  is  relatively  simple  and 
well  documented.  However,  the  design  of  step  frequency  waveforms  is  complex  due  to  the 
conflicting  effect  of  waveform  parameters  on  various  desired  performance  goals.  Thus,  in  this 
report  effort  is  made  to  systematize  the  design  process  for  the  high  PRF  step  frequency 
waveft)rm.  The  main  point  behind  the  design  of  high  PRF  waveforms  is  to  keep  the  target  in 
clutter  free  area  and  the  target  doppler  unambiguous.  This  involves  selection  of  PRF  based 
on  minimum  and  maximxun  target  velocities  and  other  parameters.  PRF  should  be  large 
enough  that  target  at  maximum  expected  velocity  stays  unambiguous  in  doppler.  The  other 
constraint  is  based  on  the  fact  that  target  doppler  from  minimum  target  velocity  exceeds  the 
clutter  doppler.  These  constrains  along  with  resolution  requirements  lead  to  the  high  PRF 
waveform  design. 

In  this  chapter,  the  method  used  to  compute  the  desired  parameters  for  the  high  PRF 
design  case  will  be  explained,  but  first  two  parameters  will  be  defined  and  the  expression  for 
constraints  on  PRF  will  be  derived. 
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1. 


Definition  of  two  constants 


A  dimensionless  parameter  P  that  relates  the  radial  target  motion  over  the  time 
duration  of  the  step  frequency  waveform  to  the  resolution  of  the  profile  can  be  defined  as: 


p  _  vNT 


(3.1) 


where: 

V  =  the  radial  velocity  between  the  radar  and  the  target  (closing  velocity); 
NT=  the  duration  of  the  step  frequency  waveform; 

AR  =  the  processed  range  resolution. 


This  is  the  ratio  of  how  fer  the  target  moved  with  respect  to  the  range  resolution  cell  during 
the  integration  time.  For  all  velocities  with  magmtude  greater  than  zero,  this  represents  a 
mismatch  in  the  DFT  (which  acts  like  a  matched  filter),  resulting  in  attenuation  and  dispersion 
of  the  HRR  profile.  In  addition  to  the  attenuation  and  dispersion  that  will  result  in  a  reduction 
in  range  resolution  and  S/N  ratio,  the  uncompensated  radial  velocity  will  shift  the  HRR 
profiles  by  L  FFT  bins,  where  L  is  given  by: 


L  = 


(3.2) 


where: 

fo  =  the  nominal  transmitter  frequency  (center  frequency); 
B  =  the  total  processed  bandwidth. 
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2.  Derivation  of  PRF  constraints 

To  get  a  moving  target  in  clutter  free  zone  after  DFT  processing,  the  motion  induced 
doppler  shift  should  be  such  that  shift  L  exceeds  the  clutter  extent, 

L  >  NxAf.  (3.3) 


From  Equations  3. 1,  3.2  and  3.3 


2NTv 

X 


k  NxAf. 


From  above  v.*,  is  obtained  as 


A  c 
2NAR 
A  cx  1 

If'm'AR 


(3.4) 


(3.5) 


The  constraint  for  maximum  velocity  can  be  derived  by  considering  the  worst  case 
scenario  for  a  target  located  at  the  end  of  the  range  gate  (of  width  ct/2)  at  the  maximum 
speed.  To  keep  the  doppler  unambiguous,  doppler  shift  and  spread  together  should  not 
exceed  the  PRF  and  therefore 

L+P  <  N-NxAf.  (3.6) 

Substituting  the  value  of  L  and  P  into  this  last  equation. 
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^NTv  vNT  ^ 

A  AR 

vT  X  ^  k  A 

- <  — -tA/ —  . 

AR  2T  2T  2T 


V  + 


(3.7) 


Last  term  in  right  hand  side  of  above  equation  can  be  recognized  as  the  minimum  velocity 
from  Equation  3.5  and  therefore 


V 


< 


A 

2T 


-  V  . 
min 


/ 

V 

max 

V 


+  V  . 
mm 


< 


(3.8) 


Now.  solving  for 


1+- 


2AR 


+  V 


(3.9) 


Equality  will  give  the  minimum  PRF  If  the  &ctor  ,i/((2AR)  «  1,  then  the  previous 

equation  can  be  reduced  to 


(3.10) 


3.  Design  Method 

The  initial  spedfications  are  the  nominal  carrier  frequency  (/,),  range  resolution  {AR), 
maximum  radial  velocity  (v„^X  minimum  radial  velocity  (v„^„ ),  time  on  target  {tot)  and  due 
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to  proces^g  constraints,  the  maximum  number  of  pulses  The  waveform  parameters 
that  we  want  to  calculate  are  the  PRF  (/.),  number  of  pulses  (N),  frequency  step  size  (Af) 
and  the  pulse  width  (r). 

The  proposed  method  consists  of  five  steps,  as  follows: 

1.  Selection  of  PRF  from  given  equation .  This  equation  gives  the  minimum  PRF  and 
the  exact  value  can  be  determined  using  an  iterative  procedure  later. 

2.  The  number  of  pulses  is  determined  from  tot  and  /.  : 


N  ^  {tot)f^ 


0 


B 


00 


fr> 


where  Og  is  the  antenna  beamwidth  and  4  ts  the  scan  rate. 

3.  The  frequency  step  size  is  calculated  from 


(3.11) 


A/  = 


c 

INAR’ 


(3.12) 


where  c  is  the  speed  of  fight  (3  x  10*  m/sec). 

4.  The  pulse  width  is  chosen  from  Equation  3.5,  which  can  be  rewritten  as 


T  <  _  - 

■  ^  fr¥' 


(3.13) 


N  and  r  will  impact  the  S/N  ratio  which  determines  the  radar  performance,  that  is 
probability  of  detection  and  probability  of  false  alarm  for  given  radar  parameters.  Therefore 
it  should  be  ascertained  from  the  following  equation  that  chosen  and  r  yield  adequate  S/N 
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ratio  to  satisfy  the  radar  performance  requirements. 


S  PjG^A^o(TiV) 
N  (4nfRHT^FL’ 


(3.14) 


where 

P,  =  the  transmitted  peak  power  (watts); 

G  =  the  antenna  gain  of  the  transmitter/receiver; 

A  =  the  wavelength  (meters); 

a=  the  radar  cross  section  of  the  target  (square  meters); 
r=  the  pulse  width  (meters); 

the  number  of  pulses  coherently  integrated  within  one  scan; 

R  =  the  detection  range  of  the  target  (meters); 
k  =  Boltzmaim’s  constant  (1.38  x  10'^  watt-second/^K); 

To  =  the  noise  temperature  (°K  ); 

F = the  receiver  noise  figure; 

Z,  =  a  loss  factor  incorporating  all  system  losses. 

By  convention.  To  is  taken  to  be  290°K,  which  is  close  to  room  temperature  (300 ‘^K) 
and  conveniently  makes  the  product  kTo  a  round  number  (4  x  10'^^  watt-second/°K). 
5.  Flexibility  in  PRF  constraint  can  be  used  to  satisfy  other  design  objectives  such  as 
the  search  range; 
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where  R„  is  the  unambiguous  range  and  Bip  is  the  intermediate  frequency  bandwidth. 
N,f„  rand  Af  are  changed  iteratively  such  that  the  required  search  range  is  obtained. 
Extent  o£R„  should  not  increase  clutter  as  range  resolution  is  already  fixed.  However 
the  number  of  ’s  will  increase  clutter  (which  should  not  matter  if  cancellation  is 
adequate). 

The  initial  spedfications  and  the  final  calculated  design  parameters  can  be  summarized 
in  the  following  table: 


Specifications  Calculated  Design  Parameters 

Range  resolution  (AR)  Minimum  PRF  (fj 

Nominal  carrier  frequency  ( f^)  Number  of  pulses  (N) 

Time-on-target  (tot)  Frequency  step  size  (Af ) 

Minimum  radial  velocity  Pulse  width 

Maximum  radial  velocity  (v^) 

Maximum  number  of  pulses  (N„^) 

Table  1.  Summary  of  speciScations  and  design  parameters. 
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4.  Graphical  implementation  of  the  Design  Method 
Using  the  equations  mentioned  in  the  previous  section,  a  graphical  method  to  design 
the  waveform  was  developed.  The  first  graph  is  a  plot  of  the  range  resolution  versus  minimum 
PRF,  for  four  different  nominal  carrier  fi'equendes. 


Figure  5 .  Minimum  PRF  versus  range  resolution  for  four  different  nominal  carrier 
frequencies. 

This  graph  was  created  using  Equation  3.9  and  it  can  be  observed  that  as  the  value 
of  the  desired  range  resolution  decreases  (corresponding  to  an  improvement  in  the  range 
resolution),  the  minimum  PRF  increases.  This  graph  is  used  to  determine  the  PRF  for  a  given 
range  resolution  and  carrier  frequency. 
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The  next  graph.  Figure  6,  is  a  plot  of  the  PRF  versus  the  number  of  pulses,  and  was 


generated  using  Equation  3.11. 


Figure  6.  Number  of  pulses  versus  PRF  for  three  different  times-on-target  (5, 10, 
and  15  ms). 


This  plot  is  used  to  determine  the  number  of  pulses  for  a  given  time-on-target  (tot) 
and  PRF,  computed  in  the  previous  step  from  Figure  5.  Then,  and  using  Equation  3.12,  the 
frequency  step  size  was  plotted  versus  the  number  of  pulses,  for  three  different  values  of 
range  resolution  (0.3,  0.5,  and  1  meter)  as  can  be  observed  in  Figure  7.  This  graph  is  used 
to  determine  the  frequency  step  for  a  given  range  resolution  and  number  of  pulses  computed 
in  the  previous  stqp  from  Figure  6.  The  frequency  step  size  required  for  the  same  number  of 
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pulses  increases  as  the  desired  range  resolution  decreases.  Therefore,  the  better  range 
resolution  is  desired,  the  higher  the  frequency  step  size  must  be. 


Figure  7.  Frequency  step  versus  number  of  pulses  for  three  different  range  resolutions 
(0.3, 0.5,  and  1  meter). 


Finally,  the  last  step  is  to  compute  the  pulse  width  using  Equation  3.13,  and  all  the  calculated 
design  parameters  from  Table  1  will  be  obtained. 

Note  that  in  all  the  calculations  and  plots  above  described  was  considered  to  be 
150  m/s  and  ,  1000  m/s.  The  value  considered  for  was  2048. 

Consider  an  example  to  determine  the  design  parameters  for  a  step  frequency  radar 
with  a  nominal  carrier  frequency  of  10  GHz,  time-on-target  of  15  ms  and  a  desired  range 
resolution  of  0.3  m.  From  Figure  5  a  minimum  PRF  of  80  KHz  would  be  obtamed.  Using  this 
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value  and  Figure  6,  N  would  be  approximately  equal  to  1200.  Now  with  N=1200,  dr=0.3m 
and  u^g  Figure  7,  the  frequency  step  would  be  approximately  0.42  MHz.  The  last  step  is  to 
compute  the  value  of  the  pulse  width  using  Equation  3.14  and  the  result  would  came  up  to 
be  t=298  ns. 

Interpolation  should  be  used  in  this  plots  whenever  the  specifications  are  different 
from  the  values  in  the  plots.  However  some  times  the  desired  values  are  beyond  the  plotted 
values.  Additionally,  not  only  we  want  to  compute  the  desired  parameters  that  satisfy  the 
initial  specifications,  but  also  we  may  want  to  satisfy  some  parameter  constraints.  Therefore 
a  MATLAB  code  was  written  in  order  to  fulfill  this  requirement.  This  program  can  be  found 
in  Appendbc  B. 


5.  Computer  implementation  of  the  Design  Method 

The  waveform  design  parameters  are  constrained  to  the  foUowing  values; 

-  Pulse  width  -  between  50  nsec  and  5000  nsec 

-  Frequency  step  -  should  be  0.4,  0.6, 0.8  or  1  MHz 

-  PRF  -  should  be  an  integer  multiple  of  5  KHz 

The  program  is  written  in  a  way  so  that  it  wiU  ask  for  the  desired  range  resolution,  nominal 
carrier  frequency,  and  the  time-on-target  and  these  data  should  be  entered  from  the  keyboard. 
With  this  input  data  the  program  wiU  compute  the  minimum  values  using  the  equations 
previously  presented.  These  values  are  iterated  until  the  constraints  mentioned  above  are  all 
satisfied  and  the  desired  design  parameters  obtained.  The  minimum  values  are  shown  in 
brackets.  The  program  will  also  compute  the  actual  range  resolution  and  time-on-target 
corresponding  to  the  rounded  values  of  the  PRF,  frequency  step,  pulse  width  and  number  of 
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pulses.  Figure  8  illustrates  how  the  screen  will  appear. 


Enter  the  desired  range  resolution  in  meters  :  0.3 

Enter  the  nominal  carrier  frequency  in  GHz  :  10 

Enter  the  time-on-target  in  msec  :  15 

The  calculated  design  parameters  are  : 

(minimum  values) 

-  PRF  =  140.00  KHz 

-  Frequency  step  =  0.40  MHz 

-  Pulse  width  =  180.00  nsec 
•  Number  of  pulses  =  2048 

-  Actual  range  resolution  =  0.19  m 

-  Actual  time-on-target  =  14.63  msec 

(80.00  KHz) 

(0.42  MHz) 

(300.00  nsec) 

(1200) 

Figure  8.  Example  of  the  screen  when  running  progrcanl.m 
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IV.  STEP  FREQUENCY  WAVEFORM  ANALYSIS  USING  THE 
AMBIGUITY  FUNCTION 

A.  INTRODUCTION 

In  this  chapter,  the  ambiguity  function  for  a  specific  step  frequency  waveform  will  be 
determined.  The  waveform  with  the  following  parameters  is  of  practical  interest: 

tj  =  0.1  psec, 

T  =  5  psec, 

Af=lMHz, 

N  =  500, 

where  t^  is  the  pulse  width,  T  is  the  PRI,  Af  is  the  frequency  step,  and  N  is  the  number  of 
pulses.  However,  it  is  difficult  to  compute  and  plot  this  ambiguity  function  with  the  desired 
resolution  because  of  the  very  large  amount  of  computations  involved.  Not  only  the  amount 
of  computations  is  large  but  also  the  plotting  is  hard  to  handle  because  the  files  can  easily 
reach  several  megabytes.  Therefore,  waveforms  with  parameters  similar  to  the  waveform  of 
interest  but  easier  to  compute  will  be  investigated,  to  bring  out  the  key  characteristics  of  the 
waveform  under  investigation.  First,  theoretical  dimensions  will  be  verified  with  a  few  cases. 
Then,  these  theoretical  dimensions  will  be  used  for  particular  case  of  interest  which  is  difficult 
to  compute  directly. 

Since  the  ambiguity  function  of  the  step  frequency  waveform  contains  elements  of 
ambiguity  functions  of  linear  frequency  modulated  (LFM)  pulses  and  train  of  constant 
frequency  pulses,  the  ambiguity  functions  of  these  waveforms  will  be  introduced  along  with 
the  general  introduction  to  ambiguity  function  in  Appendix  A. 
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B.  AMBIGUITY  FUNCTION  OF  THE  STEP  FREQUENCY  WAVEFORM 

The  transmitted  signal  of  the  step  frequency  radar  (as  shown  in  Figure  1)  can  be 
represented  mathematically  as  follows: 

s/it)  =  A^f^s(t-nT)  ,  (4,1) 

K=0 


with  s(t)  defined  as  s(t-nT)  =  1,  if  nT<  nT+t^ , 

=  0,  elsewhere,  (4.2) 

and  «  =  0, ,  N-l.  is  the  amplitude  of  the  transmitted  signal.  Equation  4. 1  can  be  rewritten 
as  follows: 


Sit) 


N-\ 

n=0 


g/2it/or 


=  A.  S(t)  e 


(4.3) 


The  ejqjression  in  brackets  is  the  complex  envelope  of  the  pulse  sequence  represented  by  S(t). 
The  received  signal  can  be  written  as 

S,(0  -  A,  S((-u)  ,  (4,4) 


and  the  ambiguity  function  as 


=  I  / 


N-\ 

E 

/M=0 


N-\ 


Y^s(t-mT)  ^  s*(t-nT-z)  e  /2’'»Ar(«-c) 


n=Q 


(45) 
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Changing  variables  ( t-mT  =  the  final  expression  for  the  ambiguity  function  of  the  step 
fi’equency  waveform  is  obtained: 

v-i  v-i 

5^  gj2%m'^Lfr  g-j2iunntifr  gj2v.nt4^z 

m-0  n-0 

(4.6) 

X  fs(t')  s‘(l'-(n-m)T-x)  dt'  \  . 

-OO 

The  theoretical  dimensions  of  the  contours  of  the  ambiguity  surface  for  the  step 
fi-equency  waveform  are  shown  in  Figure  9.  The  overall  dimensions  are  of  length  21SIT  along 
the  delay  axis  and  2NA/ along  the  firequency  axis.  It  can  be  observed  that  the  distance 
between  component  contours  along  the  delay  axis  is  equal  to  the  PRI  of  the  waveform  and 
along  the  fi-equency  axis  is  equal  to  the  inverse  of  the  PRI.  The  ambiguity  diagram  of  the  step 
fi-equency  waveform  can  be  obtained  fi-om  the  ambiguity  diagram  of  the  constant  frequency 
pulse  train  (same  parameters),  by  rotating  the  delay  axis  by  LflT.  This  can  be  stated  more 
accurately  by  saying  that  the  ambiguity  function  of  the  step  firequency  radar  is  x(t, 
where  xix,  is  the  ambiguity  function  of  the  constant  firequency  pulse  train.  Apart  fi-om 
rotation  of  horizontal  axis  (and  all  cuts  parallel  to  it)  the  individual  spikes  also  rotate.  This  is 
similar  to  linear  fi-equency  modulation  pulses.  The  extent  of  the  central  column  of  spikes  along 
the  fi-equency  axis  for  zero  delay  is  2/ts ,  value  which  decreases  when  we  move  away  fi-om  the 
center.  The  central  spike  as  shown  in  Figure  9  is  inchned  A/TT  with  respect  to  the  horizontal 
axis,  indicating  the  range-doppler  coupling  of  the  step  fi^equency  waveform.  The  projection 
of  the  central  peak  on  the  delay  axis  is  2ts ,  which  is  the  overall  delay  uncertainty.  However 
if  the  doppler  is  known,  the  delay  uncertainty  is  given  by  2/(NA^  as  indicated  in  the  figure. 
The  projection  of  the  central  peak  on  the  frequency  axis  is  also  indicated  in  the  figure.  If  the 
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Figure  9.  Level  contour  of  the  ambiguity  surface  of  a  step  frequency  waveform. 
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delay  is  known,  the  frequency  uncertainty  is  2/NT  which  is  the  same  as  for  the  constant 
frequency  pulse  train.  However  if  target  delay  is  not  known,  frequency  uncertainty  is  much 
larger  due  to  range-doppler  coupling. 

C.  AMBIGUITY  FUNCTION  OF  STEP  FREQUENCY  WAVEFORMS 

Four  waveforms  will  be  studied  in  this  section.  These  waveforms  are  defined  by  the 
following  parameters: 

1.  N=5,  t,  =1  ps,  PRI  =5  ps,  Af =1  MHz. 

2.  N=10,  t^  =1  ps,  PRI  =5  ps,  Af =1  MHz. 

3.  N=10,  f,  =0.1  ps,  PRI  =5  ps,  Af =1  MHz. 

4.  N=500,  t^  =0.1  ps,  PRI  =5  ps,  Af =1  MHz. 

1.  First  waveform 

The  results  for  this  waveform  are  shown  in  Figures  10  through  14.  The  ambiguity 
diagram  is  represented  in  Figure  10.  It  can  be  observed  that  it  is  spiky  like  the  one  for  the 
constant  frequency  pulse  train,  but  it  is  tilted  at  an  angle.  The  global  view  of  the  contour  plot 
of  the  ambiguity  function  for  this  waveform  is  shown  in  Figure  1 1(a).  The  rotation  of  the  plot 
with  respect  to  the  axes  is  obvious  in  the  figure  and  the  rotation  angle  is  given  by  Aj/7T.  The 
delay  and  frequency  axis  dimensions  match  the  theoretical  values  as  from  -NT  to  NT  and  from 
-NA/  to  NA/.  Figure  1 1(b)  is  a  close  up  view  of  the  contour  plot.  Figure  12(a)  is  a  further 
magnification  of  the  contour  plot  and  Figure  12(b)  shows  the  details  of  the  central  peak.  From 
the  contour  plot  it  is  clear  that  the  central  peak  is  a  inclined  ridge  with  the  highest  magnitude 
in  the  center  of  the  contour.  Figure  13(a)  is  a  cut  of  the  ambiguity  surface  along  the  delay 
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frequency  [MHz] 


delay  [microsec] 


Figure  10.  Ambiguity  diagram  of  a  step  frequency  waveform. 

(N=5,  pulse  width=l  microsec,  PRI=5  microsec,  Af=l  MEiz) 


frequency  [MHz] 


delay  [mIcrosec] 

(a) 


delay  [microsec] 

(b) 

Figure  1 1 .  Contour  plots  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=5,  pulse  width=l  microsec,  PRI=5  microsec,  Af=l  MHz) 

(a)  Global  view. 

(b)  Magnified  view. 


frequency  [MHz] 


delay  [microsec] 

(b) 

Figure  12.  Contour  plots  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=5,  pulse  width  =1  microsec,  PRI=5microsec,  Af=l  MHz) 

(a)  Magnified  view. 

(b)  Magnified  view  of  the  central  peak. 
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Figure  13.  Time  profiles  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=5,  pulse  width  =1  microsec,  PRI=5  microsec,  Afi=lMHz) 

(a)  Global  view. 

(b)  Magnified  view. 
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(b) 

Figure  14.  Frequency  profile  of  the  ambiguity  diagram  of  a  step  fi'equency  waveform. 
(N=5,  pulse  width  =1  microsec,  PRI=5  microsec,  Af=lMHz) 

(a)  Global  view. 

(b)  Magnified  view. 
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domain  at^=0,  which  portrays  only  one  peak.  This  suggests  that  there  is  no  range  ambiguity 
for  the  waveform  if  the  target  doppler  is  known.  Figure  13(b)  is  a  magnification  of  Figure 
13(a)  which  gives  nuU-to-null  width  of  0.4  which  is  the  same  as  computed  from  the  formula 
2/(NA/’).  Note  that  this  is  also  the  output  of  the  matched  filter  in  the  time  domain.  The  range 
resolution  of  this  waveform  is  0.2  as  given  by  1/(NA/).  Figure  14(a)  is  a  cut  of  the  ambiguity 
surface  along  the  fi-equency  domain  at  i:=0.  Theory  suggests  that  there  should  be  (2T/tj,)-l 
peaks  in  the  frequency  domain,  which  is  matched  by  the  actual  value  of  9.  Figure  14(b)  is  a 
detail  of  the  central  peak  and  should  have  a  null-to-null  width  of  2/(NT)  which  matches  with 
the  actual  value  of  0.08. 

2.  Second  waveform 

For  the  second  case,  the  waveform  has  the  same  parameters  as  in  the  previous  one 
except  that  the  number  of  pulses  is  increased  to  10.  This  would  affect  the  range  resolution 
(1/NAf)  and  the  doppler  resolution  (1/NT)  which  are  now  0.1  and  0.02  respectively. 
However,  some  other  quantities  of  interest  such  as  the  component  contours  in  the  fi’equency 
axis  stay  the  same.  The  results  can  be  observed  in  Figures  15  through  18.  Figure  15(a)  gives 
the  global  view  of  the  contour  plot  of  the  ambiguity  diagram  for  this  waveform.  Figures 
15(b),  16(a)  and  16(b)  are  successively  increasing  magnifications  of  the  contour  plot.  Figure 
16(b)  is  a  detailed  plot  of  the  central  peak  giving  a  range  resolution  and  fi-equency  resolution 
of  0.1  and  0.02,  respectively.  This  can  be  confirmed  in  Figures  17(b)  and  18(b)  which  are 
magnified  views  of  the  time  and  fi-equency  profiles,  respectively. 
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Figure  15.  Contour  plots  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=10,  pulse  width=l  microsec,  PRI=5  microsec,  A^l  MHz) 

(a)  Global  view. 

(b)  Magnified  view. 
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Figure  16.  Contour  plots  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=10,  pulse  width=l  microsec,  PRI=5  microsec,  Af=l  MHz) 

(a)  Magnified  view. 

(b)  Magnified  view  of  the  central  peak. 
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Figure  17.  Time  profiles  (cut  at  zero  fi'equency)  of  the  ambiguity  diagram  of  a  step 
fi'equency  waveform.  (N=10,  pulse  width=l  microsec,  PRI=5  microsec, 
Af=lMHz) 

(a)  Global  view. 

(b)  Magnified  view. 
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Figure  18.  Frequency  profiles  of  the  ambiguity  diagram  of  a  step  fi'equency  waveform. 
(N=10,  pulse  width  =1  microsec,  PRI=5  microsec,  A^lMHz) 

(a)  Global  view. 

(b)  Magnified  view. 
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3. 


Third  waveform 


For  this  waveform  the  pulse  width  is  reduced  to  0. 1  ps  and  the  rest  of  the  parameters 
are  the  same  as  in  the  previous  case.  The  results  can  be  observed  in  Figures  19  through  24. 
Figure  21  represents  the  central  peak  of  the  ambiguity  diagram  of  this  step  frequency 
waveform  which  is  similar  to  the  one  in  Figure  A.7. 

Now  that  all  the  theoretical  dimensions  previously  defined  were  verified,  they  can  be 
used  to  calculate  the  dimensions  of  any  particular  case  of  interest.  This  will  be  done  for  the 
fourth  and  last  case. 
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Figure  19.  Contour  plot  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=10,  pulse  width=0. 1  microsec,  PRI=5  microsec,  Af=l  MHz) 
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Figure  20.  Time  profile  (cut  at  zero  frequency)  of  the  ambiguity  diagram  of  a  step 
fi-equency  waveform.  (N=10,  pulse  width=0. 1  microsec,  PRI=5  microsec 
Af=lMHz) 
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Figure  21 .  3D  plot  of  the  central  peak  of  the  ambiguity  diagram  of  a  step  frequency 
waveform.  (N=10,  pulse  width=0.1  microsec,  PRI=5  microsec, Af=l  MHz) 
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Figure  22.  Contour  plot  of  the  central  peak  of  the  ambiguity  diagram  of  a  step 
frequency  waveform.  (N-10,  pulse  width=0.1  microsec,  PRI=5  microsec, 
AM  MHz) 
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Figure  23 .  Time  profile  (cut  at  zero  delay)  of  the  central  peak  of  the  step  fi'equency 
waveform.  (N=10,  pulse  width  =  0.1  microsec,  PRI  =  5  microsec, 
Af=lMHz) 
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Figure  24.  Frequency  profile  (cut  at  zero  delay)  of  the  central  peak  of  a  step 
frequency  waveform.  (N=10,  pulse  width=0. 1  microsec,  PRI=5  microsec, 
Af=lMHz) 
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4. 


Fourth  waveform 


Using  the  results  from  Section  B,  which  have  been  confirmed  for  the  waveforms 
discussed  so  far,  ambiguity  diagram  figures  were  sketched  for  a  waveform  with  500  pulses. 
The  amount  of  computations  required  to  generate  the  ambiguity  surface  by  computer  is 
enormous.  The  3D  arnbiguity  diagram  of  this  waveform  is  several  magnitudes  more  complex 
than  3D  pictures  of  midtown  Manhattan.  One  can  realize  the  problem  of  resolution  of  this 
figure  and  its  printing  on  a  normal  size  paper  or  even  in  a  much  larger  size  one. 

A  contour  plot  of  the  ambiguity  diagram  is  sketched  in  Figure  25.  Each  short  slant  line 
is  a  spike.  In  doppler  cBmension  the  spikes  are  PRF  (200  KHz)  apart  and  spread  over  (20 
MHZ).  In  delay  dimension  (at  fj=0)  there  are  about  ten  significant  spikes  (in  reality  there  are 
more  but  others  will  be  small  in  magnitude).  This  means  that  the  IF  filter  will  pass  ten 
frequency  lines,  which  implies  that  the  return  from  ten  different  ranges  due  to  separate  pulses 
may  arrive  at  the  same  time.  Thus,  one  may  say  that  the  target  will  be  ambiguous  in  range  and 
that  it  can  lie  in  ai^  of  the  ten  range  zones,  i^art  from  the  fact  that  the  target  range  ambiguity 
needs  to  be  resolved  it  will  also  compete  against  clutter  from  ten  different  range  zones. 

If  the  waveform  PRF  is  high,  the  doppler  from  a  target  moving  at  the  highest  expected 
velocity  may  be  kept  unambiguous.  In  such  cases,  spikes  along  the  doppler  dimension  may 
be  ignored,  otherwise  they  have  to  be  considered  for  clutter  calculations.  Figures  26  and  27 
represent  cuts  along  the  delay  and  frequency  dimensions  indicating  the  potential  range  and 
doppler  ambiguities,  respectively.  Figure  28  gives  the  detail  of  the  central  peak  (note  that 
other  peaks  or  spikes  are  of  similar  shape  but  may  be  of  different  magnitude  and  located  at 
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different  places  in  range-doppler  map).  At  =0  the  peak  is  4  ns  wide  (nuU-to-null)  in  delay 
dimension.  Similarly,  the  width  of  the  peak  in  frequency  dimension  (at  t=0)  is  given  by  2/NT 
which  comes  out  to  be  800  Hz. 
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Figure  25.  Contour  plot  of  the  ambiguity  diagram  of  a  step  frequency  waveform. 
(N=500,  pulse  width=0.1  microsec,  PRI=5  microsec,  A^l  MHz) 
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Figure  26. 


Time  profile  (cut  at  zero  fi'equency)  of  the  ambiguity  diagram  of  a  step  fi'equency 
waveform.  (N=500,  pulse  width=0.1  microsec,  PRI=5  microsec,  Af=l  MHz) 


51 


Figure  27.  Frequency  profile  (cut  at  zero  delay)  of  the  ambiguity  diagram  of  a  step  fi-equency 
waveform.  (N=500,  pulse  width=0.1  microsec,  PRI=5  microsec,  A:^l  MHz) 
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Figure  28.  Dimensions  of  the  central  peak  for  the  case  N=500,  pulse  width  =  0.1  ps, 
PRI  =  5  ps,  Af  =  1  MHz.  (Values  indicated  are  null-to-hull;  generally  3dB 
values  are  given  which  are  a  httle  less  than  half  of  the  indicated.) 
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V.  CONCLUSIONS 


Step  frequency  has  several  unique  features  and  advantages  over  conventional  high 
resolution  waveforms.  It  simultaneously  allows  pulse  compression  and  doppler  processing. 
It  can  be  implemented  on  existing  radars  with  the  addition  of  a  frequency  synthesizer  and 
suitable  signal  proces^g.  This  waveform  achieves  high  resolution  with  narrow  instantaneous 
bandwidth  and  thus  reduces  the  requirements  on  A/D  sampling  rate.  Further,  more  multiple 
time  around  clutter  is  kept  limited  by  judicious  choice  of  waveform  parameters  which  reduces 
the  amount  of  clutter  entering  the  IF  amplifier.  This  in  turn  reduces  the  requirement  for 
(fynamic  range  of  A/D.  It  is  difficult  to  achieve  high  sampling  rate  in  A/D  when  combined  with 
large  dynamic  range.  Lowering  the  requirement  on  sampling  rate  and  d3mamic  range  makes 
the  implementation  of  high  resolution  systems  feasible. 

In  this  report,  emphasis  is  put  on  the  design  of  a  waveform  for  surface  based  radars 
to  detect  small  moving  targets.  Because  of  the  many  tradeoffs  involved  in  the  design  of  a  step 
frequency  waveform,  the  choice  of  the  waveform  parameters  is  not  straight  forward.  To 
systematize  the  waveform  parameter  choice,  a  design  procedure  was  developed  for  step 
frequency  waveform.  The  proposed  method  determines  the  waveform  parameters  for  given 
radar  specifications.  Two  different  implementations  of  this  method  were  presented:  the 
graphical  and  the  computer  implementation.  The  graphical  implementation  is  a  simple  and 
quick  way  of  calculating  the  minimum  design  parameters.  The  computer  implementation  takes 
into  account  parameter  constraints  imposed  by  the  hardware. 

The  ambiguity  diagram  for  a  specific  waveform  of  interest  to  NaRD  was  developed. 
Since  the  waveform  consists  of  a  large  number  of  pulses  it  is  almost  impossible  to  directly 
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compute  the  ambiguity  fimctioa  Therefore,  waveforms  with  parameter  close  to  the  waveform 
of  practical  interest  but  easier  to  compute  were  investigated,  in  order  to  bring  out  the  key 
characteristics  of  the  waveform  of  interest.  The  theoretical  dimensions  of  the  step  frequency 
waveform  were  defined  and  verified  with  a  few  cases.  Once  this  was  done,  the  particular  case 
of  interest  was  investigated.  It  contains  the  elements  of  ambiguity  function  of  a  single  LFM 
pulse  and  a  train  of  constant  frequency  pulses.  The  3D  plot  of  the  ambiguity  function  of  a  step 
frequency  waveform  is  spiky  like  the  one  of  a  constant  frequency  pulse  train  and  is  tilted  at 
an  angle  like  the  ambiguity  diagram  of  a  LFM  pulse.  The  width  of  the  central  spike  along  the 
delay  axis  at  zero  frequency  is  2/(NA/)  as  compared  to  2t,  for  the  constant  frequency  pulse 
train,  thus  verifying  the  high  range  resolution  potential  of  this  waveforms.  The  range 
resolution  can  be  increased  by  increasing  the  product  NA/"  (and  without  decreasing  the  pulse 
width).  For  the  constant  frequency  pulse  train,  high  range  resolution  can  only  be  improved 
by  decreasing  the  pulse  width.  On  the  frequency  axis,  the  spike  width  at  zero  delay  is  equal 
to  2/(NT)  as  also  is  the  case  for  the  constant  frequency  pulse  train.  Therefore,  the  frequency 
resolution  for  the  step  frequency  waveform  is  the  same  as  for  the  constant  frequency 
waveform. 

Further  work  on  the  comparative  analysis  of  step  frequency  waveform  should  be 
performed  to  compare  it  with  conventional  high  PRF  and  medium  PRF  waveforms  for 
detection  of  small  targets.  Its  use  in  inverse  synthetic  aperture  radar  (ISAR)  for  moving  target 
identification  should  also  be  investigated. 
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APPENDIX  A 


A.  DEFINITION  OF  AMBIGUITY  FUNCTION  AND  ITS  PROPERTIES 
1.  Deflnition  of  Ambiguity  Function 

A  radar  waveform’s  ambiguity  function  is  probably  the  most  complete  statement  of 
the  waveform’s  inherent  performance.  It  is  a  formula  that  quantitatively  describes  the 
interference  caused  by  a  point  target  return  located  at  a  different  range  and  velocity  from  a 
reference  target  of  interest.  It  reveals  the  range-doppler  position  of  ambiguous  responses  and 
defines  the  range  and  doppler  resolution.  [Ref.  2,  pp.  74] 

This  quantitative  description  of  the  ability  of  a  waveform  to  resolve  two  or  more  radar 
reflectors  at  arbitrarily  different  ranges  and  velocities,  constitutes  an  important  feature  for 
quick  assessment  of  the  interference  level  with  which  a  target  of  interest  must  compete  when 
it  is  in  the  vicinity  of  other  radar  reflectors.  Although  it  is  seldom  used  as  a  basis  for  practical 
radar  tystem  deagn,  it  provides  an  indication  of  the  limitations  and  utility  of  particular  classes 
of  radar  waveforms,  and  gives  the  radar  designer  general  guidelines  for  the  selection  of 
suitable  waveforms  for  various  applications. 

The  ambiguity  fimction  of  the  waveform  s(t)  can  be  defined  in  terms  of  the  cross¬ 
correlation  of  a  doppler-shifted  version  of  the  waveform,  that  is  s(t)  exp(j27tfjt)  with  the 
unshifted  waveform.  Using  the  definition  of  cross-correlation,  it  follows  that 

-t-c» 

X  Xt-X)]  dt ,  (A.  1) 

where  t  is  the  delay  time  and_^  is  the  doppler  frequency  shift.  Rearranging  the  terms  in  the 
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integral  produces  a  common  form  of  the  ambiguity  function  as  shown  in  Equation  A.2. 


Ix(t^,4)I  =  1 1  s(t)s\t-x)  . 

—00 

A  normalized  expression  is  obtained  by  requiring  that 

+00 

fWfdi=\  . 


(A.2) 


(A3) 


W^th  this  normalization,  the  magnitude  of  the  ambiguity  function  has  a  value  of  unity  at  the 
origin. 

2.  Properties  of  the  Ambiguity  Function 

The  ambiguity  function  has  the  following  properties: 


Peak  value  of  \yi(xj^)\  =  |x(0,0)|  =  E , 

(A4) 

IX(-V/rf)|  =  \xMa)\ , 

(A.5) 

+0O 

lx  ("c,  0)1  =  1^5(0 , 

-oo 

(A6) 

+0O 

IX(0,/rf)|  =  1  j s\t)  , 

—00 

(A7) 

+00  +00 

/  j\li'^,f^\dxdf^  =  E  . 

(A8) 

Equation  A4,  states  that  the  peak  value  of  the  ambiguity  function  occurs  at  the  origin  and  it 
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is  equal  to  E,  the  energy  contained  in  the  echo  signal.  Equation  A.  5  shows  the  ambiguity 
fimction’s  symmetry.  Equation  A6  indicates  that  the  ambiguity  function  along  the  time  delay 
axis  is  the  auto-correlation  function  of  the  complex  envelope  of  the  transmitted  signal. 
Equation  A7  states  that  along  the  frequency  shift  axis  the  ambiguity  function  is  proportional 
to  the  spectrum  of  s^  (t).  Finally,  Equation  A.8  states  that  the  total  volume  under  the 
ambiguity  surface  is  a  constant  equal  to  E. 

3.  Ambiguity  Diagram 

It  is  common  to  refer  to  as  the  ambiguity  surface  of  the  waveform.  The 

shape  of  this  ambiguity  surfece  depends  entirely  on  the  waveform  parameters.  The  plot  of  the 
ambiguity  surface  is  called  ambiguity  diagram.  The  ideal  ambiguity  diagram  consists  of  a 
single  spike  of  infinitesimal  thickness  at  the  origin  and  is  zero  everywhere  else.  The  single 
center  spike  eliminates  any  ambiguities,  and  its  infinitesimal  thickness  at  the  origin  permits  the 
frequency  and  the  echo  delay  time  to  be  determined  simultaneously  to  an  high  degree  of 
accuracy.  It  also  pennits  the  resolution  of  two  targets  no  matter  how  close  together  they  are 
on  the  ambiguity  diagram.  Naturally,  this  ideal  diagram  does  not  exist.  The  two  reasons  for 
this  can  be  found  in  the  properties;  first,  and  accordingly  with  Equation  A.4,  the  mavimnni 
hdght  of  the  ambiguity  function  is  E  and  secondly  the  volume  under  the  surface  must  be  finite 
and  equal  to  E,  as  stated  by  Equation  A.8.  However,  a  reasonable  approximation  is  given  in 
Figure  Al.  This  ambiguity  function  only  has  one  peak  and  therefore  does  not  cause  any 
ambiguity.  However  the  single  peak  might  not  be  narrow  enough  to  satisfy  the  requirements 
of  accuracy  and  resolution.  If  the  single  central  peak  is  made  too  narrow,  it  may  cause  other 
smaller  peaks  to  occur  in  regions  other  than  the  origin,  and  therefore  cause  ambiguities.  The 
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Figure  A.  1.  An  approximation  to  the  ideal  ambiguity  diagram  “After  Ref.  [1].” 

requirements  for  accuracy  and  unambiguity  are  not  always  possible  to  satisfy  simultaneously. 

The  particular  waveform  transmitted  by  a  radar  is  chosen  to  satisfy  the  requirements 
for  (1)  detection,  (2)  measurement  accuracy,  (3)  resolution,  (4)  ambiguity,  and  (5)  clutter 
rejection.  The  ambiguity  diagram  may  be  used  to  assess  how  weU  a  waveform  can  achieve 
these  requirements.  Each  of  these  will  be  discussed  briefly. 

The  requirements  for  detection  do  not  place  any  demands  on  the  shape  of  the 
transmitted  waveform  except  that  it  be  possible  to  achieve  with  practical  radar  transmitters, 
and  the  maximum  value  of  the  ambiguity  function  is  an  indication  of  the  detection  capabilities 
of  the  radar.  The  accuracy  with  which  the  range  and  the  velocity  can  be  measured  by  a 
particular  waveform  depends  on  the  width  of  the  central  spike  along  the  time  and  frequency 
axis.  The  resolution  is  also  related  to  the  width  of  the  central  spike,  but  in  order  to  resolve 
two  closely  spaced  targets  the  central  spike  must  be  isolated.  It  cannot  have  any  high  peaks 
nearby  that  can  mask  another  target  close  to  the  desired  target.  A  waveform  that  yields  good 


resolution  will  also  yield  good  accuracy,  but  the  reverse  is  not  always  so. 

A  continuous  waveform  (a  single  pulse)  produces  an  ambiguity  diagram  with  a  single 
peak.  A  discontinuous  waveform  can  result  in  peaks  in  the  ambiguity  diagram  at  other  values 

The  pulse  train  is  a  common  example.  The  presence  of  additional  spikes  can  lead  to 
ambiguity  in  the  measurement  of  target  parameters.  An  ambiguous  measurement  is  one  in 
which  there  are  several  choices  available  for  the  correct  value  of  a  parameter,  but  only  one 
choice  is  appropriate.  Thus  the  correct  value  is  uncertain.  The  ambiguity  diagram  permits  a 
visual  indication  of  the  ambiguities  possible  with  a  particular  waveform.  The  ambiguity 
problem,  detection  and  accuracy  are  related  to  a  single  target,  whereas  resolution  applies  to 
multiple  targets. 

The  ambiguity  diagram  may  be  used  to  determine  the  ability  of  a  waveform  to  rgect 
clutter  by  superimposing  on  the  t,  plane  the  regions  where  clutter  is  found.  If  the 
transmitted  waveform  is  to  have  good  clutter-rejection  properties  the  ambiguity  function 
should  have  little  or  no  response  in  the  regions  of  clutter. 

The  problem  of  S3mthesizing  optimum  waveforms  based  on  a  desired  ambiguity 
diagram  specified  by  operational  requirements  is  not  normally  feasible.  The  approach  to 
iselecting  a  waveform  with  a  suitable  ambiguity  diagram  is  not  systematic  but  rather  by  trial 
and  error. 

The  name  ambiguity  function  is  somewhat  misleading  since  this  function  describes 
more  about  the  waveform  than  just  its  ambiguity  properties.  This  name  was  given  to  this 
function  in  order  to  demonstrate  that  the  total  volume  under  it,  is  a  constant  equal  to  E, 
independent  of  the  shape  of  the  transmitted  waveform.  Thus  the  total  area  of  ambiguity,  or 
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uncertainty,  is  the  same  no  matter  how  the  ambiguity  surface  is  distributed  over  the 
plane.  [Ref.  1,  pp.  418-420] 


B.  AMBIGUITY  FUNCTION  OF  A  SINGLE  PULSE 

The  unmodulated  single  pulse  is  widely  used  in  old  generation  radars  for  search  and 
track  functions  (as  magnetrons  are  not  easily  modulated)  and  where  range  accuracy  and 
resolution  requirements  can  be  met  with  a  pulse  wide  enough  to  provide  sufficient  energy  for 
detection.  It  has  the  minimum  ratio  of  time  sidelobe  extent  to  compressed  pulse  width  and  is 
used  in  inexpensive  radars  where  signal  generation  and  processing  costs  must  be  minimized. 
The  single  pulse  of  a  sine  wave  can  be  defined  as 

,  (A.9) 


where  s(t)  is  the  complex  envelope  of  the  signal,  defined  as  follows: 

s(t)  =  1,  if  0<  t<  ts, 

=  0,  elsewhere,  (A.  10) 

with  ts  equal  to  the  pulse  width. 

Using  Equation  A.2,  the  ambiguity  fimction  of  the  single  pulse  can  be  written  as 
follows: 


(i_m1 

sin[nr^(l-|t|//^)/J 

i  “J 

zero  elsewhere.  (A.  1 1 ) 
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The  ambiguity  diagram  of  the  single  pulse  is  shown  in  Figure  A.  2  and  the  contour 
plot  can  be  observed  in  Figure  A.3.  Profiles  of  the  ambiguity  function,  taken  through  its  peak, 
are  shown  in  Figures  A.4  and  A.  5.  The  time  profile  is  a  triangle  (the  autocorrelation  of  a 
rectangle  is  a  triangle)  with  a  half-voltage  width  equal  to  the  pulse  width.  The  frequency 
profile  is  a  sine  function  whose  main  lobe  width  is  twice  the  reciprocal  of  the  pulse  width. 
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Figure  A.  3. 
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Contour  plot  of  the  ambiguity  diagram  of  a  single  pulse, 
(pulse  \vidth=l  microsec) 
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Figure  A.4  Time  profile  of  the  ambiguity  diagram  of  a  single  pulse, 
(pulse  width  =1  microsec) 
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Figure  A.  5.  Frequency  profile  of  the  ambiguity  diagram  of  a  single  pulse. 

(pulse  width  =1  microsec) 
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C.  AMBIGUITY  FUNCTION  OF  THE  LEVEAH  FREQUENCY  MODULATED 
PULSE 

The  linear  frequency  modulated  (LFM)  pulse  is  commonly  used  to  increase  range 
accuraty  and  resolution  when  long  pulses  are  required  to  get  reasonable  signal-to-noise  ratios 
(10  to  20  dB).  This  waveform  can  be  used  for  detection  of  targets  with  unknown  velocity 
since  the  doppler  sensitivity  is  low. 

The  LFM  pulse  can  be  represented  mathematically  as 


S(t)  =  s(t)e 


(A.12) 


where  s(t)  is  the  same  as  defined  in  Equation  A.  10,  and  k  is  the  rate  of  frequency  change  in 
Hz/sec.  The  ambiguity  function  of  the  LFM  pulse  can  be  written  as  in  Equation  A.  13. 


+00 


(i-Wl 

sin[7tr^(l-|T|/^^)(/^+AT)] 

I  d 

if  I'cl  ^  ( 


zero  elsewhere.  (A  13) 


Comparing  this  equation  with  Equation  A.2,  we  can  see  that  they  are  identical,  except  that 
j^is  replaced  by  fj-  kt.  Therefisre ,  we  can  conclude  that  the  ambiguity  function  of  the  LFM 
pulse  is  a  shifted  version  of  the  ambiguity  fimction  of  a  single  pqlse.  This  relation  is 
represented  graphically  in  Figure  A.  6  and  it  can  be  seen  that  the  ambiguity  function  of  the 
LFM  pulse  is  just  a  rotated  version  of  the  ambiguity  function  of  the  single  pulse. 


67 


fi 


Figure  A.6.  Contour  comparison  for  the  single  pulse  and  the  LFM  pulse. 


Figures  AJ,  A8,  A9,  and  A 10,  represent  the  ambiguity  diagram,  contour  plot,  and 
time  and  frequency  profiles  (taken  through  the  peak),  respectively. 
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Figure  A.  7.  Ambiguity  diagram  of  a  LFM  pulse, 
(pulse  width=l  microsec) 
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Figure  A.  9.  Time  profile  of  the  ambiguity  diagram  of  a  LFM  pulse, 
(pulse  width  =1  microsec) 


Figure  A.  10.  Frequency  profile  of  the  ambiguity  diagram  of  a  LFM  pulse. 

(pulse  width  =1  microsec) 


71 


D.  AMBIGUITY  FUNCTION  OF  THE  CONSTANT  FREQUENCY  PULSE  TRAIN 

Practical  radars  employ  waveforms  consisting  of  constant  frequency  pulse  trains  and 
therefore  it  is  important  to  study  their  ambiguity  functions.  The  ambiguity  functions  for  these 
pulse  trains  reveal  ambiguous  responses  in  range  and  doppler.  Performance  for  specific 
surveillance  applications  can  be  understood  in  terms  of  unambiguous  range-doppler  regions 
of  operation  determined  by  radar  pulse  repetition  frequency,  pulse  duration,  and  pulse 
bandwidth.  For  high  resolution  jq)plications,  the  ambiguity  surface  of  individual  pulses  of  the 
pulse  train  is  also  of  interest.  [Ref  2,  pp.  74-75] 

For  single  pulse  delay  and  frequency  measurement  accuracies  depend  on  the  single 
parameter  of  pulse  width.  With  the  pulse  train  this  situation  can  be  avoided.  The  delay 
accuracy  depends  on  the  pulse  width  as  before,  but  the  frequency  accuracy  is  now  determined 
by  the  total  duration  of  the  pulse  train.  Therefore,  both  accuracies  are  independent  of  one 
another.  The  price  that  has  to  be  paid  for  this  capacity  of  independently  controlling  delay  and 
frequency  accuracies  is  that  additional  peaks  occur  in  the  diagram,  which  in  turn  cause  range 
and  doppler  ambiguities.  In  practice,  the  radar  designer  tries  to  select  the  PRI  in  order  to 
make  all  targets  of  interest  appear  in  the  vicinity  of  the  central  peak,  and  all  the  other  peaks 
occur  as  far  from  this  region  as  possible. 

Most  radars  use  this  t5q)e  of  waveform,  which  can  be  mathematically  represented  as 
follows: 


S(t)  =  ,  (A.  14) 

w=0 


where  is  the  complex  envelope  of  the  single  pulse  of  the  transmitted  signal,  defined  as 
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s(t-nT)  =  \,  if  nT<  t<nT+t^, 

=  0,  elsewhere.  (A.  15) 


In  this  equation,  Nvs,  the  number  of  pulses,  T  is  the  PRI  and  4  is  the  pulse  width.  The 
ambiguity  function  of  the  constant  frequency  pulse  train  can  be  written  as 

lx('C,/rf)l  =  If  dt\  .  (A.  16) 

^  W=0  772=0 


Changing  variables  ( t-nT=t  ’ ),  Equation  A.  16  becomes 

v-i  v-i  . 

ix('c,/rf)l  =  lE  E  fs(t')  s\t'-{m-n)T-x)  dt'\  .  (a.17) 

72=0  772=0  ^ 

-00 

Figure  A 1 1  is  a  level  contour  of  the  ambiguity  surface  of  a  constant  frequency  pulse 
train.  This  figure  shows  that  the  width  of  the  central  peak  along  the  delay  axis  is 
approximately  twice  the  pulse  width,  and  that  along  the  frequency  axis  is  2/NT.  The 
interpeak  distance  along  the  delay  axis  is  equal  to  the  PRI,  and  along  the  frequency  axis 
1/PRI,  which  is  the  PRF  of  the  waveform.  The  total  non-zero  extent  along  the  delay  axis  is 
2NT  and  along  the  frequency  axis  is  Figure  A 12  represents  an  example  of  an  ambiguity 

diagram  for  the  constant  frequency  pulse  train.  The  contour  plots  in  Figures  A 13  and  A 14 
show  that  the  results  obtained  agree  with  the  theoretical  ones  previously  defined.  This  can 
also  be  observed  in  Figures  A 15  and  A 16,  which  represent  the  time  and  frequency  profiles, 
respectively.  The  time  profile  of  the  central  pulse  is  a  triangle  and  the  frequency  profile  is  a 
sine  function,  as  expected. 
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Figure  A12.  Ambiguity  diagram  of  a  constant  frequency  pulse  train. 

(N=5,  pulse  width=l  microsec,  PRI=5  niicrosec,  AF=0) 
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Figure  A.  13.  Contour  plots  of  the  ambiguity  diagram  of  a  constant  frequency  pulse 
train  (N=5,  pulse  width=l  microsec,  PRI=5  microsec,  Af=0) 

(a)  Global  view. 

(b)  Magnified  view. 


Figure  A.  14.  Contour  plots  of  the  ambiguity  diagram  of  a  constant  frequency  pulse 
train.  (N=5,  pulse  width  =1  microsec,  PRI=5  microsec,  AfN)) 

(a)  Magnified  view. 

(b)  Magnified  view  of  the  central  peak. 
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Figure  A.  1 5.  Time  profiles  of  the  ambiguity  diagram  of  a  constant  fi'equency  pulse 
train.  (N=5,  pulse  width  =1  microsec,  PRI=5  microsec,  Af=0) 

(a)  Global  view. 

(b)  Magnified  view. 
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Figure  A.  16.  Frequency  profile  of  the  ambiguity  diagram  of  a  constant  frequency 
pulse  train .  (N=5,  pulse  width  =1  microsec,  PRI=5  microsec,  Af=0) 

(a)  Global  view. 

(b)  Magnified  view. 
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APPENDIX  B.  MATLAB  SOURCE  CODES 


%  File  name;  programl.m 
%  Title:  Plot  generation 
%  Date  of  last  revision;  4  Mar  1 996 
%  Written  by:  Paulo  A.  Soares 

%  Comments;  This  program  generates  3  plots  for  the  SFWF  design. 


clear 

% - specifications - 

vmin=150;vmax=1000; 

f[)l=3e9;fD2=5e9;fD3=7e9;fD4=10e9; 

totl=5e-3;tot2=10e-3;tot3=15e-3;tot4=20e-3; 


% - computations  &  plots - 

c=3e8; 

lambdal=c^  1  ;lambda2=c/ro2; 
lambda3=c/fD3  ;lambda4=c^04; 

.% - plot  #1 - 

di=0.01:0.01:l; 

fi’l=(dr.\'vmax)+((2/lamb(lal)*(vmax+vmin)); 

fi:2=(dr.\vmax)+((2/lambda2)*(vmax+vmin)); 

fi‘3=(dr.\vmax)+((2/lambda3)*(vmax+vmin)); 

fi-4=(dr.\vmax)+((2/lambda4)*(vmax+vmin)); 

figure(l) 

plot(dr^  l/le3, 'w')4iold  on 

plot(drJ&2/le3/w')4iold  on 

plot(dr3/le3,’w')Jiold  on 

plot(dr^4/le3,V) 

axis([0.01,l,10,110]) 

xlabelCRange  Resolution  (meters)') 

ylabelCMinimum  PRF  (Khz)') 

text(0.75^7,'fo=3GHz'),text(0.75,42,'fo=5GHz') 

text(0.75,57,'fo=7GHz'),text(0.75,81,'fo=10GHz') 

text(0.2,100,'fo = nominal  carrier  fiequency') 

lambda=0.03 :0.00 1 :0. 129; 
k=(dr.\vmax); 

firl=(lambda.\2)*(vmax+vmin)+k; 

Nll=totl*fi-l; 

N21=tot2*fi;l; 

N31=toG*fi-l; 

figure(2) 

plot(fi:l/le3,Nl  l,V)4iold  on 
plot(frl/le3,N21,'wOJiold  on 
plot(fi-l/le3,N3  l,'w')Jtiold  on 
axis([20,180,0,3000]) 
xlabel(TRF  (KHz)') 
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ylabel(TSr) 

text(140,900/tot=5ins'),text(140,1700,'tot=l(hiis') 
text(140^550,tot=l  Sms') 
text(50,2500,’tot  =  time-on-target') 

o/o - .plot  #3 - 

deltarl=0.3; 

deltar2=0.5; 

deltar3=l; 

N=l:  1:2048; 

deltafl=(2*N*deltarl  *  le6).\c; 
deltaf2=(2*N*deltar2*  le6).\c; 
deltaf3=(2*N*deltar3  *  le6).\c; 
figiire(3) 

plot(N,deltafl,'w')Jbold  on 
plot(N,deltaf2,’w’)4iold  on 
plot(N,deItaf3,’w’) 

axis([0^100Al]) 

xlabel('N') 

ylabel(Trequency  Step  (MHz)') 

text(1700,0.12,'dr=lm'),t«£t(1700,0.2i;dr=0.5m'),text(1700,0.33;dr=0.3m') 
text(1200,0.8,'dr  =  range  resolution') 
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%  Filename;  program2.m 
%  Title:  Main  program  for  SFWF  design 
%  Date  of  last  revision:  28  Feb  1996 
%  Written  by:  Paulo  A.  Soares 

%  Comments:  This  program  computes  the  exact  design  parameters  fimin,  df  N  and  tan,  and  then 
%  rounds  it  in  accordance  with  previously  defuied  criteria. 

clear 

c=3e8;Nmax=2048; 

vmin=150jvmax=1000; 

dr=input('  Enter  the  desired  range  resolution  in  meters :  ■); 
while  dr  <=  0 

dr=input(’  You  must  enter  a  positive  value.  Please  try  again :  •) 
end 

fD=input(’  Enter  the  nominal  carrier  frequency  in  GHz :  *); 
while  fO  <=  0 

fO=input(’  You  must  enter  a  positive  value.  Please  try  again : ') 
end 

lambda==c/(fr)*  1  e9); 

tot=input('  Enter  the  time-on-target  in  msec : '); 
while  tot  <=  0 

tot=input('  You  must  enter  a  positive  value.  Please  try  again :  ’) 
end 

frinin=(((2/lambda)*(vmax+vmin))+(vmax/dr))/le3; 
fr=ceil(frmin); 
while  rem(fr,5)  ~  0 
fi=fr+l; 
end 

Nmin=ceil(tot*fr); 

N=Nmin; 

ifN<=2 

N=2; 

end 

elseif  N=3  |  N=4 
N=4; 
end 

elseif  N>4  &  N<^ 

N=8; 

end 

elseif  N>9  &  N<^17 
N=16; 
end 

elseif  N>17  &  N<=34 
N=32; 
end 

elseif  N>34  &  N<=68 
N=64; 
end 

elseif  N>68  &  N<=136 
N=128; 
end 

elseif  N>136  &  N<=272 
N=256; 
end 
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elseifN>272&N<=544 

N-512; 

end 

elseif  N>544  &  N<=1088 
N=1024; 
end 

elseif  N>1088 
N=Nmax; 
end 

Nroimd=N; 
fiTOiind=NroiindAot; 
fiTOund=ceil(fiTound); 
while  rem(fiTOiind,5) 
fiToimd===firromid+l ; 
end 

totroimd=Nromid/&Tonnd; 

dfinin=c/(2*Nmin*dr*  le6); 

df=c/(2*Nroiind*dr*  le6); 

stepsizedf=0.2; 

dfroimd=ceil(d£^stepsi2edf); 

dJ&oiind==dfix)und*stepsizedf; 

taiimin=(2*\inin)/(lambda*dfQiin*fiimn); 

tau=(2*vmin)/(lambda*djfroimd*fiTomid); 

stepsizetau=10; 

taiiroimd=ceil(tau/stepsizetau); 

tauronnd=tauround*stepsizetau; 

acdr=(c/(2*Nround*dfround))*  le-6; 

stepsizeacdr=0.0 1 ; 

acdiToimd===ceil(acdr/stepsizeacdr); 

acdiToimd=acdrTomid*stepsizeacdr; 

j|^rin1f('\n  The  calculated  design  parameters  are:W) 

:fyrintfC  (minimum  values  in  brackets)\nW) 

^rintfC  -PRF  =  %3.2fKH2  (%3.2fKHz)\n’^ound,fimin); 

^rintfC  -Frequency  step  =  %3.2f MHz  (%3.2f MQHLz)W,dfix)und,dfrnm); 

^rintfC  -Puke  width  =  %3.2f  nsec  (%3.2fnsec)\n^tauround,taumin); 

fyrintfC  -Number  of  pulses  =  %4.0f  (%4.0f)WJ^round,Nmin); 

^rintf(’  -Actual  range  resolution  =  %3.2f  m\n',acdrround); 
fyrintf(*  -Actual  time-on-target  =  %3 .2f msec  ',totround); 


%  Filename:  prograin3*iii 

%  Title:  Ambiguity  function  of  the  simple  pulse  by  correlation  method 
%  Date  of  last  revision:  6  Apr  1996 
%  Written  by:  Paulo  A.  Soares 

%  Comments:  This  program  draws  four  plots  of  the  ambiguity  function  for  a  single  pulse. 


clear; 

tau=l; 

nx=50; 

ny=201; 

t=[0:nx-l]*tau/nx; 

*=t(2)-t(l); 

z=D; 

fd=liiispace(- 1 0/tau,  1 0/tau^y); 

ul=ones(l4ix); 

form=laiy 

u2=ul  *exp(j*2*pi*fd(m)*t); 
c=xcorr(u2,ul).*dt; 
z=[z;abs(c)]; 
end 

Hfliplr(-t),t(2:nx)]; 
figure(l) 
mesh(t/d^),grid 
xlabel('delay  [microsec]') 
ylabel(’fiequency  [NMz]') 
zlabel('magnitude') 

£i«ure(2) 

contour(t/d^,7),grid 
x]abel('delay  [microsec]') 
ylabel(’frequency  [MHz]') 
figure(3) 

plot(tX(iiy+l)/2,;)>') 
xlabel('delay  [microsec]') 
ylabelCmagnitude') 
figure(4) 

plot(fd^(:^),'w’) 
xlabel('£requency  [MHz]') 
ylabel('magnitude') 
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%  Filename:  program4.m 

%  Title:  Ambiguity  function  of  the  LFM  pulse  by  correlation  method 
%  Date  of  last  revision:  6  Apr  1996 
%  Written  by:  Paulo  A.  Soares 

%  Comments:  This  program  draws  four  plots  of  the  ambiguity  function  for  a  LFM  pulse. 


clear, 

tau=l; 

mu=4; 

nx=50; 

ny=201; 

t=[Onix-l]*tau/nx; 

dt=t(2)-t(l); 

z=D; 

fd=liiispace(-10/tau,  lO/tai^ny); 
form=l3iy 

ul=ejq)(j*pi*mu.*t.*t); 

u2=coiij(ul).*ejq)(j*2*pi*fd(m)*t); 

c=conv(ii2iliplr(ul)).*dt; 

z=[z;abs(c)]; 

end 

t=[fliplr(.t),t(2:nx)]; 
figure(l) 
mesh(^d;z),grid 
xlabel('delay  [micFOsec]') 
ylabelCfiequency  [MHz]') 
zlabelCmagnitude’) 
figure(2) 

contour(t^d^5),grid 
xlabelCdelay  [microsec]') 
ylabel('firequency  [MHz]') 
figure(3) 

plot(t;z((ny+l)^,:),'w') 
xlabel('del^  [microsec]') 
ylabel('magnitude') 
figure(4) 

plot(fd^(:^);w') 
xlabelCfiequency  [MHz]') 
ylabelCmagnitude') 
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%  Filename;  program5.m 

%  Title:  Ambiguity  limction  of  the  step  frequency  radar  by  the  correlation  method 
%  Date  of  last  revision:  26  Apr  1996 
%  Written  by:  Paulo  A.  Soares 

%  Comments:  This  program  computes  and  plots  several  figures  of  the  ambiguity  fimction  of  the 
%  step  fi^quency  radar,  using  the  correlation  method. 


O/o- - 

% 

%- - 

Parameters  initialization 

clear^ 

%  clears  all  variables 

tau=l; 

%  tau=pulse  widfii 

T=5; 

%  T=period  or  pulse  repetition  interval  (PRI) 

N=10; 

%  N=niimber  of  pulses 

dfi=l; 

%  dP=frequency  step 

sampfreq=40; 

%  sampfreq=sampling  frequency 

NT=N*T; 

%  NT=total  period  of  the  waveform 

iisamptau=sampfreq*tau; 

%  nsamptau=number  of  samples  in  one  pulse  width 

nsamptO=sampfreq*(T-tau); 

%  nsamptO=number  of  samples  in  the  interpulse  period 

nsamptau=floor(nsamptau); 

nsamplO=floor(iisamptO); 

iisampT=iisamptau-httsainptO; 

%  nsampT=number  of  samples  in  one  PRI 

%. - 

% 

%. - 

Definition  of  axis  and  resolution 

P=[ones(  1  ,nsamptau),zeFos(  1  ^isamptO)] ; 

t=iiiispace(0,NT,N*nsampT); 

%  time  axis  definition 

ny=301; 

%  number  of  points  in  the  freq.  axis 

fd=lmspace(-5/T,5/T^y); 

%  freq.  axis  definition 

dt=t(2)-t(l); 

%  dt=time  interval 

% 

Computations 

0/^ - 

s=D; 

m=l:l:N; 

forn=l:N 

s=[s,ejq)(j*2*pi*(m(n)-l)*df*t(nsampT*(n-l)+l:nsampT*n)).*P]; 

end 

z=D; 

forL=l:ny 

Md(L); 

sl=conj(s).*exp(j*2*pi*f*t); 

s2=s; 

c=conv(s  1  Thplr(s2)).  *dt; 

%z=[z;(abs(c)).^2]; 

27=[z;abs(c)]; 

end 

[l,r]=size(z); 

t=[fliplr(-t),t(2:N*nsampT); 
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Plotting  of  the  results  and  axis  labeling 


O/o. - 

% 

o/o. - 

figure(l) 
niesh(^d^),giid 
xlabelCdelay  [microsec]') 
ylabel('frequency  [MHz]') 
figuie(2) 
contour(t^d^5) 
axis([-10  10  -1  1]) 
xlabelCdelay  [microsec]') 
ylabelCfreqnency  [MHz]') 
figure(3) 

plot(t^((ny+l)/2,:)) 
axis([-30  30  0  12]) 
xlabel('del^  [microsec]') 
ylabel('magnitude') 
figure(4) 

plot(fd'^(:,(r)-l)/2),'w') 
xlabel('freqnency  [MHz]') 
ylabel('magnitude') 
figure(5) 
contour(t^d^,5) 
axis([-1.5  1.5  -0.2  0.2]) 
xlabelCdelay  [microsec]') 
ylabelCfiequency  [MHz]') 
figure(6) 

plot(t^((ny+l)/2,:)) 
axis([.1.5  1.5  0  12]) 
xlabelCdelay  [miCTosec]') 
ylabelCmagnitude') 
figure(7) 

plot(fd'^(;,(rfl)/2)') 
xlabelCficquency  [MHz]') 
axis([-0.12  0.12  06]) 


%  Filename;  program6.m 

%  Title:  Ambiguity^  function  of  the  step  frequency  radar  by  the  equation  method 
%  Date  of  last  revision:  26  Apr  1996 
%  Written  by:  Paulo  A.  Soares 

%  Comments:  This  program  computes  and  plots  several  figures  of  the  ambiguity  function  of  the 
%  step  frequency  radar,  using  the  equation  me&od. 

%- - 

%  Parameters  initialization 

O/o- - 

clear; 
tau=l; 

T=5; 

N=10; 
df=l; 
d=tau/T; 


o/o. - 

% 

Definition  of  axis  and  resolution 

nx=101; 

%  nx=nnmber  of  points  in  the  time  axis 

t=linspace(-N*T,N*Tjnx); 

%  time  axis  definition 

ny=3001; 

%  ny=number  of  points  in  the  frcq.  axis 

fd=lmspace(-5/T,5/T  ^y); 

%  fi'eq.axis  definition 

% 

Computations 

o/o- -  - 

p=floor(t/T); 

r=(t/T)-p; 

z=[]; 

form=l;ny 

x=n; 

f=fd(m); 

forn=l:nx 

pp=p(n); 

rr=r(n); 

fl=f-pp*df; 

C=f-(2*pp+rr)*df; 

f3=f-(pp+l)*df; 

f4=f-(2*pp+n+l)*df; 

ifrr>=0&rr<d 

iffl=0 

xl=l; 

else 

X 1  =(sin(pi*fl  *(tau-rr*T)))7(pi*fl  *(tau-rr*T)); 
end 

if  rem(f2*T,l)==0 
x2=N-abs(pp); 
else 

x2=(sin(pi*f2*(N-abs(pp))*T))7(sin(pi*f2*T)); 

end 

je=[x,(tau-rr*T)  *xl  *x2]; 
elseif  rKl  &  rr>=(l-d) 
iff3==0 
xl=l; 


%  clears  all  variables 
%  tau=pulse  width 
%  T=pulse  repetition  interval  (PRI) 
%  N=number  of  pulses 
%  df=fiiequency  step 
%  d=duly  cycle 
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else 

xl=(siii(pi*G*(tau-Krr-l)*T)))./(pi*f3*(tau+(rr-l)*T)); 

end 

if  rem(f4*T,  1)=0 
x2=N-abs(pp+l); 
else 

x2=(sm(pi*f4*(N-abs(pp+l))*T))./(sin(pi*f4*T)); 

end 

x=[x,(tau+(rr-l)*T)*xl*x2]; 

else 

x=[x,0]; 

end 

end 

z=[z;x]; 

end 

z=abs(z); 

O/o- - 

%  Plotting  of  results  and  axis  labeling 

o/^ - 

figure(l) 
mesh(t;fd;z),grid 
xlabel('delay  [microsec]') 
ylabel('j&equency  [MHz]') 
figure(2) 
contour(t^<l^,5) 
axis([-6  6-0.2  0.2]) 
xlabel('delay  [microsec]') 
ylabel('fiequen(gf  [MHz]') 
figure(3) 

plot(U((ny+l)/2,:)) 
xlabel('delay  [microsec]') 
ylabel('magnitiide') 
figure(4) 

plot(fd'^(:,(nx+l)/2)') 
xlabel(^fi'equency  [MHz]') 
ylabel('magnitude') 
figure(5) 
contoiir(t/d^,5') 
axis([-0.7  0.7  -0.05  0.05]) 
xlabel('delay  [mi<aosec]') 
ylabelCfiequency  [MHz]') 
figiire(6) 

plot(t^((iiy-i-l)/2,:)) 
axis([.1.5  1.5  0  12]) 
xlabel('delay  [microsec]') 
ylabel('magnitude') 
figure(7) 

plot(fd';z(:,(nx+l)/2)') 
axis([-0.05  0.05  0  0.5]) 
xlabel('fiequency  [MHz]') 
ylabel('magnitude') 
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